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EVALUATION  OF  AFTERBURNING  CESSATTO^ MECHANISMS  IN  FUEL  RICH 
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Air  Force  Research  Laboratory 
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Abstract 

A  computational  study  was  conducted  to  identify  fun¬ 
damental  physical  processes  governing  the  cessation  or 
shutdown  of  the  afterburning  of  fuel  rich  rocket  exhaust 
with  the  atmosphere.  Several  mechanisms  that  have 
been  proposed  were  examined  which  are:  1)  a  relami- 
narization  phenomenon,  2)  a  Damkohler  number  effect 
and  3)  a  classical  flame  extinction  mechanism.  Analy¬ 
sis  of  the  simulation  results  revealed  the  relaminanza- 
tion  mechanism  to  be  implausible  while  the  Damkohler 
number  effect  and  the  flame  extinction  mechanisms  were 
found  to  be  valid.  The  extinction  mechanism  was  also 
found  to  dramatically  alter  the  emission^  characteristics  ✓ 
and  enhance  the  shutdown  behavior  which  has  impor¬ 
tant  implications  with  respect  to  radiative  heat  trans¬ 
fer  to  the  body  and  missile  intercept  systems.  This  is 
a  significant  finding  because  strain  rate  induced  extinc¬ 
tion  is  a  previously  unrecognized  phenomena  occurring 
in  rocket  exhaust  plumes  during  afterburning  cessation. 


1  Introduction 

Future  development  of  booster  systems  will  be  substan¬ 
tially  influenced  by  the  temperature  environment  which 
is  encountered  in  the  body’s  base  region.  Heat  trans- 
fer  to  the  body  results  from  both  convective  and  la- 
diative  phenomena.  For  fuel  rich  systems  exhibiting 
strong  afterburning  of  the  engine  exhaust  species  with 
the  atmosphere,  the  radiative  component  of  the  base 
heat  transfer  rate  may  be  the  dominant  feature  (Reijasse 
and  Delery,  1994).  As  altitude  or  time  after  launch  is  in¬ 
creased,  the  radiative  heat  transfer  rate  will  eventually 
show  a  large  drop  in  magnitude  (Kramer,  1970),  asso¬ 
ciated  with  the  shutdown  of  afterburning  in  the  plume. 
The  time  (or  altitude)  and  rate  of  decay  of  the  afterburn¬ 
ing  shutdown  event  will  determine  the  total  heat  transfer 
to  the  body  and  thus  the  design  criteria  for  the  compo¬ 
nents  in  the  base  region.  The  radiative  emissions  during 
the  shutdown  event  is  also  relevant  to  missile  defense 
systems.  Consequently,  an  accurate  characterization  of 
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afterburning  shutdown  or  cessation  events  is  required. 

Accurately  characterizing  exhaust  plume  afterburning 
and  afterburning  cessation  is  a  difficult  problem  due  to 
the  complexities  of  the  physical  phenomena  occurring  in 
these  flows.  For  example,  strong  turbulence-chemistry 
interactions  may  exist  in  the  high  speed  plume  shear 
layer  between  the  exhaust  and  external  flow.  These  in¬ 
teractions  not  only  affect  species  production  and  tem¬ 
perature  distributions  but  may  also  impact  the  radia¬ 
tion  predictions  (Pearce  and  Varma,  1981).  In  addition, 
a  complex  shock  wave  pattern  (Dash  et  al,  1980)  exists 
in  the  plume  to  equilibrate  the  exhaust  flow  pressure 
with  the  atmosphere.  Reflected  shocks  may  excite  vorti¬ 
cal  roll  up  of  the  shear  layer  (e.g.,  Norris  and  Edwards, 
1997)  enhancing  mixing  and  may  also  provide  a  temper¬ 
ature  rise  to  stabilize  afterburning.  These  difficulties  are 
also  exacerbated  by  3-D  effects  in  multiple  nozzle  sys¬ 
tems. 

Additional  complications  arise  in  exhaust  plume 
analysis  due  to  uncertainties  with  regard  to  the  combus¬ 
tion  and  mixing  processes  occurring  within  the  engine 
itself.  First,  the  injector  pattern  and  any  acoustic  baf¬ 
fles  within  the  engine  set  up  fuel-oxidizer  striations  which 
may  not  be  completely  homogenized  before  exiting  into 
the  plume  due  to  incomplete  mixing.  Second,  fuel  film 
cooling  and  regenerative  cooling  of  the  fuel  and  oxidizer 
may  have  a  significant  impact  on  the  species  concentra¬ 
tion  and  temperature  distributions.  Third,  the  spectral 
distribution  of  the  turbulent  kinetic  energy  within  the 
engine  nozzle  has  also  been  shown  to  strongly  effect  af¬ 
terburning  and  shutdown  characteristics  (Schadow  et  al., 
1990).  Fourth,  the  presence  of  particulates  may  have  a 
significant  impact  on  the  plume  thermal  structure  and 
flowfield  (Lee  and  Netzer,  1994;  Dash,  1993)  which  in 
turn  will  affect  afterburning  and  shutdown.  Deficiencies 
in  characterizing  these  uncertainties  will  inhibit  accurate 
predictions  of  plume  afterburning  and  shutdown,  result¬ 
ing  in  poor  predictions  of  radiative  emissions  and  heat 
transfer  to  the  body. 

The  character  of  the  afterburning  shutdown  events  has 
been  observed  to  vary  among  different  missile  systems. 
For  example,  during  shutdown  some  systems  exhibit  a 
gradual  drop  in  total  radiant  intensity  over  a  wide  alti- 
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tude  range  while  other  systems  shutdown  very  rapidly 
over  a  much  narrower  range.  Application  of  the  engi¬ 
neering  level  modeling  techniques  to  predict  afterburning 
cessation  has  been  relatively  successful  for  systems  which 
exhibit  the  gradual  drop  off  type  of  shutdown  event. 
However,  for  systems  showing  abrupt  shutdown,  engi¬ 
neering  models  have  been  less  successful.  To  improve  the 
accuracy  of  the  engineering  codes,  a  basic  understand¬ 
ing  of  the  physical  mechanisms  producing  afterburning 
cessation  must  be  gained.  Such  an  understanding  is  a 
prerequisite  for  future  development  efforts  to  capture 
the  essential  physics  without  unnecessarily  compromis- 
ing  computational  efficiency. 

Several  mechanisms  have  been  proposed  to  explain 
afterburning  shutdown  in  fuel  rich  exhaust  plumes 
with  altitude  increase  which  include  1)  a  relaminar- 
ization/mixing  transition  phenomenon,  2)  a  Damkoh- 
ler  number1  effect  and  3)  a  classical  flame  blow- 
off/extinction  mechanism.  The  relaminanzation  mech¬ 
anism  is  based  on  the  assertion  that  the  shear  layer 
Reynolds  number  (Re)  reduction  with  increasing  altitude 
will  cause  the  mixing  transition  (Breidenthal,  1981)  to  be 
reached  so  that  afterburning  may  not  be  sustainable  due 
to  a  reduced  molecular  mixing  rate.  A  Damkohler  num¬ 
ber  effect  producing  afterburning  cessation  results  from 
the  competition  between  turbulent  mixing  and  chemical 
heat  release  and  radical  production.  Turbulent  mixing 
continually  cools  the  hot  engine  exhaust  through  mixing 
with  the  atmosphere  while  combustion  processes  release 
heat  and  produce  radical  species  required  to  ignite  the 
flow.  Afterburning  will  not  be  sustainable  when  turbu¬ 
lent  mixing  can  cool  the  plume  at  a  higher  rate  than 
chemical  processes  can  produce  heat  and  a  radical  pool 
sufficient  to  initiate  burning.  For  the  classical  flame  ex¬ 
tinction  mechanism,  decreasing  reaction  rates  with  in¬ 
creasing  altitude  may  cause  the  plume  shear  layer  flame 
to  become  locally  extinguished  due  to  high  turbulent 
mixing  (or  strain)  at  the  small  spatial  scales.  This  will 
result  in  a  delay  in  the  onset  of  afterburning  to  a  signifi¬ 
cant  distance  downstream  of  the  engine  nozzle  exit  plane. 
This  delay  distance  will  increase  with  altitude  until  the 
shear  layer  flame  eventually  blows-off  (i.ev  afterburning 
shutdown). 

The  purpose  of  this  study  was  to  computationally  eval¬ 
uate  the  plausibility  of  the  aforementioned  mechanisms 
to  produce  afterburning  cessation  in  fuel  rich  exhaust 
plumes.  This  has  been  accomplished  through  a  paramet¬ 
ric  study  of  the  afterburning  characteristics  of  a  generic 
Amine  booster.  This  study  considered  only  a  single  noz¬ 
zle  configuration  to  assess  these  mechanisms  before  a 
more  complex  3-D  configuration  is  studied.  To  reduce 
the  parameter  space  of  the  study,  it  was  assumed  that 

lThe  Damkohler  number  is  the  ratio  of  the  plume  characteristic 
mixing  and  chemical  time  scales  and  will  be  defined  in  a  later 
section. 


no  particulates  were  present  in  the  flow  and  that  the 
previously  mentioned  ‘real’  engine  effects  (e.g.,  mixing 
inefficiency,  regenerative  cooling,  etc.)  were  negligible. 

The  organization  of  this  paper  is  as  follows.  The  next 
section  discusses  the  computation  methodology  used  to 
carry  out  the  simulations  of  the  parametric  study.  Sec¬ 
tion  3  describes  the  application  of  this  methodology  to 
the  rocket  plume  simulations  considered  here.  Section  4 
outlines  the  parametric  study  and  presents  results.  Sec¬ 
tion  5  presents  conclusions  drawn  from  this  study. 

2  Computational  Methodology 

Turbulence  near  the  rocket  body  and  in  the  plume  causes 
the  flowfield  to  be  highly  unsteady  and  contain  a  wide 
range  of  turbulence  length  scales  due  to  the  high  Reynold 
numbers  typically  encountered  in  this  type  of  flow.  To 
accurately  capture  this  unsteady  behavior,  an  unsteady 
simulation  methodology  such  as  large-eddy  simulation 
(LES)  is  required.  However,  simulating  the  entirety 
of  the  rocket  body,  base  region  and  plume,  as  con¬ 
sidered  here,  precludes  use  of  unsteady  simulation  ap¬ 
proaches  such  as  LES  due  to  grid  resolution  require¬ 
ments  and  CPU  availability.  Consequently,  a  steady 
state,  Reynolds  averaged  approach  has  been  adopted  for 

this  study.  , 

The  body  and  plume  flowfields  were  analyzed  using 

the  structured,  compressible  Navier-Stokes  (NS)  flow 
solver  “The  General  Aerodynamic  Simulation  Program” 
(GASP)  (Walters  et  al.,  1989).  This  code  solves  the  time 
averaged  NS  equations,  including  turbulence  modeling, 
and  is  sufficiently  general  to  accommodate  chemically  re¬ 
acting  flows  and  complex  physical  geometries.  The  code 
includes  standard  models  for  thermodynamic  properties 
and  detailed  models  for  molecular  transport  properties. 
However,  the  present  simulations  only  assumed  constant 
Schmidt  and  Lewis  numbers  since  molecular  transport  is 
not  resolved  in  turbulent  simulations  using  the  Reynolds 
average  approach.  GASP  achieves  good  computational 
efficiency  through  use  of  a  fully  parallelized  shared  mem¬ 
ory  implementation  of  the  flow  solver  and  also  through 
the  use  of  a  parabolized  form  of  the  Navier-Stokes  equa¬ 
tions  (PNS)  for  parabolic  regions  of  the  flow. 

One  major  drawback  of  the  GASP  code  is  the  lack  of 
any  treatment  of  the  effect  of  turbulent  fluctuations  on 
the  mean  reaction  rate  in  the  species  conservation  equa¬ 
tions.  The  mean  reaction  rate  is  calculated  based  on 
the  mean  temperature  and  species  and  neglects  all  tur¬ 
bulent  fluctuations.  Properly  accounting  for  turbulence- 
chemistry  interactions  in  combustion  problems  has  been 
the  subject  of  intense  research  for  at  least  the  last  two 
decades.  It  is  well  known  that  neglecting  turbulent  fluc¬ 
tuations  in  the  mean  reaction  rate  terms  can  have  signifi¬ 
cant  consequences  on  the  accuracy  of  a  simulation.  This 
deficiency  precludes  the  use  of  GASP  to  evaluate  the 
previously  mentioned  flame  extinction  mechanism  in  ai- 
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terburning  plumes.  However,  this  extinction  mechanism 
was  evaluated  using  a  separate  code  to  be  described  in 

b  later  section.  .  . . 

Radiation  calculations  were  decoupled  from  the  flow- 
field  simulations  and  were  carried  out  as  described 
by  Ludwig  and  Malkmus  (1981)  and  Nelson  (1987). 
Within  this  approach,  radiation  transport  equa  ions 
were  solved  using  a  band  model  for  the  gas  phase  ab¬ 
sorption/emission,  details  of  which  are  omitted  here  for 
brevity.  The  calculations  for  the  present  study  used  a 
wide  band  pass  to  encompass  emissions  from  U"  'C  U’ 
COo  and  H20.  Total  radiation  intensity  predictions 
were  also  made  using  a  field  of  view  large  enough  to 
encompass  the  entire  plume. 


3  Rocket-Plume  Simulations 

The  simulations  considered  in  this  study  were  carried 
out  for  a  generic  Amine  booster.  Figure  3.1  presents  a 
schematic  of  the  axisymmetric  geometric  configuration 
The  exhaust  nozzle  was  assumed  to  have  an  area  ratio  ot 
10  and  protrude  aft  of  the  booster  base  as  seen  in  Figure 
3.1.  The  entire  body-base-plume  was  simulated  using 
GASP’s  NS  and  PNS  solvers.  The  elliptic  regions  at  the 
booster  nose  and  in  the  base  region  were  computed  using 
the  NS  solver  while  the  flow  over  most  of  the  body  and  in 
the  plume  was  treated  using  the  PNS  option.  Each  ax¬ 
isymmetric  solution  used  ~  25,000  grid  points  to  resolve 
the  body-base  region  and  -  60, 000  -  75, 000  points  for 
the  plume,  depending  on  the  plume  extent.  The  grids 
for  each  altitude  were  manually  adapted  around  regions 
of  high  gradients  to  ensure  proper  resolution  of  the  How 

features.  ,  , 

The  inflow  boundary  conditions  at  the  nozzle  exit 

plane  were  specified  from  a  separate  nozzle  flow  calcula¬ 
tion.  This  solution  was  carried  out  using  GASP  s  PNb 
solver  with  Chien’s  (1982.)  k-e  turbulence  model.  The 
inflow  boundary  condition  for  this  simulation  was  spec¬ 
ified  from  Mach  line  conditions  at  the  throat  generated 
using  the  TDK  package  (Nickerson,  et  al,  1989).  A  nine 
species-ten  reaction  step  mechanism  for  H2/CU  oxida¬ 
tion  was  used  for  this  calculation  as  well  as  in  the  plume. 
Figures  3.2-3.4  present  the  radial  variation  of  the  major 
species,  velocity  and  temperature,  respectively,  at  the 
nozzle  exit  plane.  For  this  case  the  flow  is  fairly  uni¬ 
form  across  the  exit.  The  temperature  (Figure  3.4)  is 
~  1330  K  in  the  core  with  a  much  higher  value  at  the 
nozzle  wall  due  to  adiabatic  boundary  conditions  im¬ 
posed  there.  The  fuel  rich  conditions  of  the  engine  are 
evident  in  Figure  3.2  where  there  is  an  excess  of  CO  and 
some  H2  with  02  virtually  depleted.  This  nozzle  solu¬ 
tion  was  assumed  to  be  invariant  with  respect  to  altitude 
and  was  used  for  all  the  body-base-plume  calculations. 

The  body-base-plume  simulations  were  carried  out  by 
specifying  the  inputs  to  the  GASP  code  as  follows.  The 
convective  flux  terms  were  discretized  using  the  Van  Leer 


flux  split  scheme  with  3rd  order  MUSCL  extrapolation 
and  min-mod  and/or  ENO  limiting.  The  viscous  terms 
within  GASP  are  discretized  to  second  order  accuracy. 
Turbulence  was  modeled  using  GASP’s  implementation 
of  the  compressibility  corrected  k-e  model.  Flow  along 
the  body  was  assumed  to  be  turbulent  and  a  wall  func¬ 
tion  was  used.  The  turbulent  Schmidt  and  Prandtl  num¬ 
bers  were  specified  as'.. 7  and  ;9,  respectively. 

Thermodynamic  properties  for  the  simulations  were 
specified  from  standard  polynomial  fits  except  at  the 
highest  altitude  which  required  the  use  of  an  equilibrium 
statistical  mechanics  approach  due  to  low  temperature  in 
the  plume  expansion  region.  Also,  molecular  transport 
properties  were  specified  using  Sutherland  curve  fits  of 
viscosity  and  thermal  conductivity  for  each  species  and 
the  species  diffusion  coefficients  were  specified  assuming 
a  constant  Schmidt  number  oL.7. 


4  Results  and  Discussion 

In  order  to  assess  the  physical  mechanisms  discussed 
in  Section  1  which  may  be  responsible  for  afterburn¬ 
ing  shutdown  in  fuel  rich  rocket  plumes,  a  parametric 
study  was  conducted.  This  study  considered  parame¬ 
ters  such  as  turbulence  level  and  chemistry  effects.  The 
different  cases  considered  in  the  study  are  listed  in  Ta¬ 
ble  4  1  The  calculations  were  carried  out  using  GASF  s 
implementation  of  the  k  -  £  turbulence  model  and  the 
previously  mentioned  9  species-10  reaction  CO/H 2  oxi¬ 
dation  mechanism.  The  turbulent  mixing  for  Case  D  was 
enhanced  by  increasing  the  production  coefficient  for  the 
turbulent  kinetic  energy  in  the  k  -  e  model  by  25%. 

Table  4.1 

Body-Base-Plume  Simulations 
Case  Turbulence  Chemistry" 


A 

B 

C 

D 


yes 

yes 

no 

yes,  enhanced 


finite  rate 
frozen 
finite  rate 
finite  rate 


Experimentation  was  used  to  determine  the  altitude 
range  over  which  afterburning  shutdown  occurred  for 
this  generic  booster  system,  given  an  assumed  trajec¬ 
tory  profile.  Within  this  range,  simulations  were  con¬ 
ducted  at  altitudes  of  25,  30  and  35  km  for  each  case 
listed  in  Table  4.1.  The  Reynolds  and  Mach  numbers 
based  on  freestream  conditions  and  the  body  radius  were 
Re  =  4  x  106  ,  2  X  106  and  1  x  106  and  Moo  -  2.6,  3.2 
and” 3.9  for  the  altitudes  of  25  ,  30  and  35  km,  respec- 

tively.  f  , 

Figures  4.1  and  4.2  present  contour  plots  of  the  tem¬ 
perature  and  C02  mole  fraction,  respectively,  for  Case  A 
at  25  km.  The  features  evident  in  this  solution  are  typi¬ 
cal 'of  all  cases.  A  curved  “barrel”  shock  is  evident  in  ig- 
ure  4.1(b)  and  is  characteristic  of  highly  underexpanded 
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plume  flows  (Dash  et  al.,  1980).  This  shock  reflects  off 
the  axis  of  symmetry  and  interacts  with  the  surrounding 
shear  layer  between  the  exhaust  gases  and  the  external 
flow.  As  seen  in  Figure  4.1(b),  the  reflected  barrel  shock 
serves  to  thicken  the  shear  layer  and  enhance  turbulent 
mixing  as  well  as  increase  the  shear  layer  temperature. 
This  shock-shear  layer  interaction  enhances  burning  an 
acts  as  an  anchor  for  the  plume  flame.  This  may  be  seen 
in  Figure  4.2  which  shows  a  rapid  rise  in  C02  produc¬ 
tion  at  this  point.  This  behavior  is  consistent  with  the 
work  of  Norris  and  Edwards  (1997)  who  applied  a  much 
more  general  unsteady  LES  technique  to  simulate  high 
speed  reacting  exhaust  flows.  Further  downstream  the 
plume  is  seen  to  be  vigorously  afterburning  as  appar¬ 
ent  from  the  spatial  distribution  of  C02  in  Figure  4.2. 
These  figures  also_  shows,  that  in  the  far  field  the  plume  ■' 
flovT  closely  resembfes  what  is  expected  for  a  standard 
turbulent  jet  flame. 

With  increasing  altitude,  the  point  where  C02  pro¬ 
duction  begins  to  rise  gradually  moves  further  down¬ 
stream  for  all  cases.  This  is  illustrated  in  Figure  4.3 
which  is  a  radial  plot  of  C02  mole  fraction  for  Case  A 
at  a  fixed  streamwise  location  x/xrefl  in  the  Plume  for 
each  altitude.  The  length  scale  xrefi  is  the  distance  be¬ 
tween  the  nozzle  exit  plane  and  the  barrel  shock  reflec¬ 
tion  point  at  the  axis  of  symmetry.  This  reference  length 
was  found  to  approximately  scale  the  results  with  respect 
to  altitude  so  that  a  fixed  x/xrefi  value  corresponds  to 
the  same  relative  location  in  the  plume  for  each  case. 
As  seen  in  Figure  4.3,  the  peak  C02  value  in  the  shear 
layer  at  25  km  is  significantly  higher  than  in  the  core 
(y  —  0),  indicating  vigorous  afterburning.  With  increas¬ 
ing  altitude,  ignition  is  delayed,  resulting  in  lower  peak 
C02  at  30  km.  At  35  km  the  plume  at  this  location  is 
extinguished.  Comparison  of  Case  A  and  B  results  show 
that  at  35  km  the  plume  is  mostly  extinguished  but  is 
still  burning  to  a  small  extent  in  the  farfield. 

The  drop  off  of  C02  production  with  respect  to  alti¬ 
tude  for  the  various  cases  correspondingly  results  in  a 
drop  off  of  the  radiation  emission  as  can  be  seen  in  Fig¬ 
ure  4.4.  This  plot  presents  contour  plots  of  the  spatial 
variation  of  the  radiant  intensity  for  Cases  A,  C  and  D. 
The  left  edge  of  each  image  in  this  figure  corresponds 
to  the  nozzle  exit  plane  and  the  downstream  direction 
is  to  the  right.  The  spatial  extent  shown  in  Figure  4.4 
encompasses  the  entire  region  of  significant  radiation  foi 
all  cases  except  Case  C.  For  this  case,  the  lack  of  tur¬ 
bulent  mixing  causes  the  plume  to  be  extraordmai  lly 
long  and  extend  outside  of  the  plotted  domain  shown 
in  Figure  4.4.  Species  and  temperature  plots  for  Case  C 
(not  shown)  reveal  that  this  case  produces  vigorous  after¬ 
burning  for  each  altitude.  However,  the  point  at  which 
product  formation  begins  to  rise  rapidly  moves  further 
downstream  with  altitude  increase,  as  evident  in  Figure 
4.4.  The  other  cases,  however,  are  nearly  extinguished 


at  35  km. 

The  apparent  drop  in  radiation  intensity  seen  in  Fig¬ 
ure  4.4  can  be  quantified  in  terms  of  the  total  radiant 
intensity  which  is  the  integration  the  spatial  field  over 
the  entire  field  of  view  for  each  case.  The  total  intensity 
is  plotted  in  Figure  4.5  with  all  values  being  normalized 
by  the  result  of  Case  A  at  25  km.  From  this  figure,  Case 
A  total  intensity  is  seen  to  drop  an  order  of  magnitude 
over  this  altitude  range  and  is  close  to  the  nonreacting 
value  of  Case  B  by  35  km.  The  emissions  of  Case  B 
result  mainly  from  the  hot  core  flow  downstream  of  the 
barrel  shock  reflection  point  and  from  high  temperatures 
in  the  high  speed  shear  layer  resulting  from  viscous  en¬ 
ergy  dissipation. 

Similar  to  the  results  for  Case  A,  Case  D  with  en¬ 
hanced  turbulent  mixing  also  exhibits  an  order  of  mag¬ 
nitude  drop  in  total  intensity  over  the  altitude  range 
considered.  However,  the  Case  D  results  are  significantly 
lower  for  each  altitude  compared  with  Case  A.  This  vari¬ 
ation  of  radiation  emissions  with  respect  to  turbulent 
mixing  rate  has  been  recognized  by  Dash  et  al.  (1980) 
in  the  context  of  turbulence  model  assessment  for  plume 
flows.  The  increased  mixing  rate  in  Case  D  also  increases 
the  drop-off  rate  of  the  total  intensity  over  Case  A.  This 
indicates  an  interaction  between  the  turbulent  mixing 
rate  and  the  chemistry,  since  a  simple  axial  scale  reduc¬ 
tion  of  the  plume  extent  due  to  an  increased  turbulent 
mixing  rate  alone  would  only  produce  a  scale  factor  re¬ 
duction  in  total  intensity  in  Figure  4.5. 

For  Case  C,  the  lack  of  turbulent  mixing  at  all  the 
simulated  altitudes  yields  a  thin  shear  layer  with  high 
temperatures  due  to  viscous  energy  dissipation.  Also, 
the  centerline  temperature  downstream  of  the  shock  sys¬ 
tem  levels  out  at  approximately  900,  850  and  750  K  for 
25,  30  and  35  km  altitudes,  respectively.  Consequently, 
both  the  hot  core  gases  and  the  high  temperature  sheai 
layer  provide  an  ignition  source  for  the  plume,  resulting 
in  the  vigorous  afterburning  and  the  high  total  intensi¬ 
ties  observed  in  Figure  4.5  for  this  case.  These  ignition 
sources  will  be  present  in  the  plume  until  much  higher 
altitudes  are  reached  (required  to  cool  the  core  gases) 
and  the  shear  is  reduced  across  the  layer  (i.e.,  velocity 
matching  of  the  core  and  external  flow  along  the  flight 
trajectory).  Consequently,  a  purely  laminar  plume  will 
burn  at  much  higher  altitudes  than  a  turbulent  one  since 
molecular  mixing  alone  is  much  less  effective  than  tur¬ 
bulence  at  cooling  the  plume  through  mixing  with  the 
atmosphere. 

The  steady  state  simulations  for  Case  C,  however, 
likely  overpredict  the  temperature  increase  in  the  plume 
shear  layer  due  to  viscous  energy  dissipation  during  tran¬ 
sition  as  altitude  is  increased.  During  transition,  the 
shear  layer  will  still  be  dominated  by  large  scale  struc¬ 
tures  (Breidenthal,  1981)  which  will  thicken  the  layer 
reducing  mean  shear  below  what  will  be  realized  in  a 
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purely  laminar  case.  Thus  viscous  energy  dissipation 
would  be  reduced,  possibly  removing  one  of  the  ignition 
sources  mentioned  in  the  previous  paragraph.  However, 
small  scale  turbulent  structures  which  promote  increased 
scalar  gradients  and  enhance  molecular  mixing  are  in 
short  supply  during  transition  (Breidenthal,  1981).  Con¬ 
sequently,  as  the  plume  shear  layer  mixes  with  the  hot 
core  gases  it  would  likely  still  ignite  due  to  the  lack  of 
small  scale  turbulence  to  quickly  homogenize  tempera¬ 
ture  gradients.  The  laminarized  plume  should  therefore 
still  burn  until  an  altitude  is  reached  which  sufficiently 
cools  the  plume  core  flow  by  expansion. 

Figure  4.6  presents  a  plot  of  the  Reynolds  number 
along  the  plume  shear  layer  for  Case  A  at  35  km.  In 
this  figure  the  Reynolds  number  is  given  by, 

ReSl  ^pAUSJJi  (1) 

where  AU  is  the  velocity  difference  across  the  layer  and 
the  overbar  represents  the  average  of  the  layer  edge  con¬ 
ditions.  The  layer  thickness  6i  is  the  distance  between 
the  points  lh  -)AAU  and  U2  +*IAU  where  Di  and  U2" 
are  the  high  and  low  speed  layer  edge  velocities  respec¬ 
tively.  The  oscillations  seen  in  this  figure  downstream 
of  x/xrefi  ~  1-5  result  from  compression  and  expan¬ 
sion  wave  interactions  with  the  shear  layer  downstream 
of  the  barrel  shock  reflection  point.  From  the  com¬ 
pressible  shear  layer  experiments  of  Goeble  and  Dutton 
(1991),  the  layer  should  be  transitioned  to  a  fully  devel¬ 
oped  turbulent  flow  by  ReSl  =  1  x  105.  IVom  Figure 
4.6,  this  value  of  ReSi  is  reached  very  close  to  the  noz¬ 
zle  exit  plane.  However,  this  transition  value  of  Re&l 
does  not  account  for  pressure  gradient  and  curvature  ef¬ 
fects  present  in  plume  flows.  Even  so,  the  plume  shear 
layer  Reynolds  number  becomes  quite  high  further  down¬ 
stream  and  can  therefore  be  expected  to  be  fully  tur¬ 
bulent.  Based  on  Figure  4.6  and  the  arguments  of  the 
previous  paragraph,  the  relaminarization/ mixing  tran¬ 
sition  mechanism  to  produce  afterburning  shutdown  of 
fuel  rich  exhaust  plumes  is  implausible. 

Another  mechanism  discussed  in  Section  1  which  may 
be  responsible  for  afterburning  shutdown  is  a  Damkoh- 
ler  number  effect.  The  Damkohler  number,  Da,  is  the 
ratio  of  representative  mixing  and  chemical  time  scales 
(i.e.,  Da  =  Tmix/Tchem )  and  characterizes  the  competi- 
tion  between  plume  cooling  due  to  turbulent  mixing  and 
plume  heating  due  to  chemical  heat  release  and  com¬ 
bustion.  Afterburning  in  exhaust  plumes  will  not  be 
sustained  when  turbulent  mixing  can  cool  the  plume  at 
a  rate  higher  than  chemical  processes  can  produce  heat 
and  a  radical  pool  sufficient  to  initiate  burning.  A  char¬ 
acteristic  mixing  time,  Tmix,  may  be  calculated  as  the 
time  of  flight  of  a  particle  traveling  from  the  nozzle  to  the 
downstream  location  where  the  shear  layer  reaches  the 
axis  of  symmetry.  A  characteristic  chemical  time,  rchem  > 
may  be  calculated  in  a  similar  manner  as  used  by  Mun- 
gal  and  Frieler  (1988),  given  the  mean  nozzle  species  and 


ambient  conditions.  As  the  altitude  is  increased,  Tchem 
will  increase  due  to  falling  reactant  concentrations.  The 
mixing  time  will  also  increase  due  to  lowering  ambient 
pressure  and  Reynolds  numbei. 

Figure  4.7  presents  a  plot  of  Da  as  a  function  of  alti¬ 
tude  for  Case  A,  The  falling  value  of  Da  seen  in  this  plot 
indicates  that  rchem  increases  more  quickly  with  altitude 
than  Tmix.  This  falling  Damkohler  number  explains  why 
afterburning  shuts  down  in  the  present  simulations.  It 
also  explains  the  increased  standoff  distance  of  vigor¬ 
ous  afterburning  seen  for  the  cases  in  Figure  4.4.  For 
Case  D  with  enhanced  mixing,  Da  will  be  lower  than 
for  Case  A  resulting  in  the  increased  drop  off  rate  of  the 
total  intensity  as  seen  in  Figure  4.5.  This  Damkohler 
number  mechanism  for  afterburning  shutdown  will  gen¬ 
erally  produce  gradual  drop  off  rates  of  plume  radiation 
as  observed  in  present  simulations.  However,  the  drop 
off  rate  produced  by  this  mechanism  will  be  dependent 
on  the  ignition  characteristics  of  the  fuel  species  within 
the  plume. 

It  should  be  noted  that  this  Damkohler  number  mech¬ 
anism  not  only  produces  the  afterburning  shutdown  ob¬ 
served  in  the  present  simulations  but  it  is  also  the  only 
shutdown  mechanism  modeled  within  GASP  and  most 
commercially  available  codes.  This  results  from  the  as¬ 
sumption  that  the  effect  of  turbulent  fluctuations  on  the 
mean  reaction  rate  are  negligible  as  discussed  in  Section 
2.  Consequently,  the  flame  extinction  shutdown  mech¬ 
anism  introduced  in  Section  1  cannot  be  evaluated  di¬ 
rectly  using  GASP  because  this  phenomena  is  a  result 
of  the  interaction  of  turbulence  and  chemistry  at  small 
spatial  scales. 

To  evaluate  the  flame  extinction  mechanism  for  af¬ 
terburning  shutdown,  an  overlay  procedure  was  devel¬ 
oped  based  on  an  elementary  flame  extinction  model 
(Sloan  and  Sturgess,  1994)  applicable  to  high  speed 
flows.  This  model  was  applied  by  Sloan  and  Sturgess 
(1994)  in  the  context  of  a  low  speed  coannular  step  com¬ 
bustor  and  is  based  on  the  Eddy  Dissipation  Concept 
(EDC)  turbulent-chemistry  interaction  model  of  Mag- 
nussen  and  Hjertager  (1976).  The  EDC  is  a  distributed 
reaction  model  (or  thick  flame  model)  and  assumes  the 
fine  scale  turbulent  structures,  within  which  combustion 
occurs,  to  be  homogeneously  mixed.  Combustion  within 
these  structures  is  then  represented  by  a  constant  pres¬ 
sure,  adiabatic  Perfectly  Stirred  Reactor  (PSR)  model 
(Glarborg  et  al.,  1986).  The  residence  time,  r,  of  this 
“reactor”  is  specified  as  the  local  hydrodynamics  time 
scale  of  the  micro-scale  structures,  which  is  the  inverse 
of  the  small  scale  strain  rate,  s.  This  strain  rate  is  esti¬ 
mated  from  the  turbulence  model  as, 

where  e,  v  and  C„,  are  the  turbulent  kinetic  energy  dis¬ 
sipation  rate,  the  laminar  mixture  viscosity  and  a  cali- 
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bration  coefficient,  respectively.  Given  this  strain  rate  or 
time  scale,  the  local  fluid  mixture  is  allowed  to  react  in 
the  PSR.  If  the  local  strain  rate  is  too  high  (or  the  time 
scale  too  small),  no  combustion  will  occur  and  the  mix¬ 
ture  is  extinguished.  For  lower  strain  rates,  the  mixture 
will  ignite  and  combustion  products  will  be  formed.  This 
model  is  in  effect  a  Damkohler  number  based  mechanism 
similar  to  what  is  currently  employed  within  GASP  by 
assumption  (as  discussed  earlier).  However,  Da  for  this 
model  is  based  on  the  micro-scale  time  as  opposed  to 
GASP  which  implicitly  uses  a  large  scale  mixing  time. 

This  simplistic  extinction  model  provides  a  binary 
on/off  switch  for  the  chemical  source  term  in  the  species 
conservation  equations.  Given  the  local  turbulent  strain 
rate  and  fluid  mixture,  the  model  uses  the  PSR  calcu¬ 
lation  to  determine  if  burning  is  extinguished.  If  extin¬ 
guished.  the  mean  reaction  rate  in  the  species  conserva¬ 
tion  equation  should  be  set  to  zero,  otherwise  chemical 
reactions  are  allowed. 

Other  more  comprehensive  extinction  models  do  ex¬ 
ist  such  as  the  thin  flame  model  of  Peters  and  Williams 
(1983).  However,  models  of  this  type  are  inapplicable  in 
their  current  form  to  the  type  of  high  speed  flows  encoun¬ 
tered  here.  These  models  represent  the  local  thermo¬ 
chemical  state  by  a  collection  or  “library”  of  laminar  op¬ 
posed  flow  diffusion  flames  given  an  assumed  joint  proba¬ 
bility  distribution  function  of  mixture  fraction  and  flame 
stretch  parameter  (Sanders  and  Lamers,  1994),  coupled 
with  an  extinction  criteria.  These  models  require  apri- 
ori  knowledge  of  local  flow  conditions  to  generate  the  li¬ 
brary  of  opposed  flow  flames.  However,  local  conditions 
are  unknown  in  the  present  context  because  there  are 
large  variations  in  pressure  and  temperature  in  both  the 
fuel  and  oxidizer  streams  due  to  high  compressibility  and 
viscous  energy  dissipation  within  the  plume  shear  layer. 
Also,  nonuniformities  in  the  nozzle  exhaust  stream  pre¬ 
clude  the  definition  of  a  single  mixture  fraction  variable. 
As  a  result,  the  simplistic  EDC  based  extinction  model 
appears  to  be  one  of  the  few  such  models  applicable  to 
the  present  flow  regime. 

To  apply  the  EDC  based  extinction  model  to  high 
speed  plume  flows  would  require  integrating  the  PSR 
calculations  within  the  main  code.  However,  this  can  be 
prohibitive  due  to  the  increased  computational  expense 
of  the  PSR  calculations  as  the  extinction  limit  is  ap¬ 
proached  as  noted  by  Sloan  and  Sturgess  (1994).  There¬ 
fore,  since  the  purpose  of  this  investigation  was  to  eval¬ 
uate  the  validity  of  the  flame  extinction  mechanism  to 
produce  afterburning  shutdown  and  not  to  make  quan¬ 
titative  predictions,  the  EDC  based  extinction  model 
was  implemented  in  an  uncoupled  overlay  calculation. 
A  code  was  written  to  make  an  uncoupled  calculation  of 
the  binary  switch  value  for  the  extinction  model  given  a 
converged  flowfield  solution.  For  this  uncoupled  analy¬ 
sis  to  be  meaningful,  the  calculation  must  be  made  on 
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a  nonreacting  flowfield  because  extinguished  regions  of 
the  flow  with  high  local  strain  produce  the  nonreacting 
solution. 

Figure  4.8  presents  a  contour  plot  of  the  binary  switch 
for  Case  B  at  altitudes  of  25,  30  and  35  km.  The  stream- 
wise  axis  for  each  altitude  has  been  scaled  by  Xrefi-  For 
this  calculation,  Cs  =  0.03  as  will  be  discussed  later.  The 
black  regions  of  the  plot  for  each  altitude  represent  the 
“on”  value  for  the  chemistry  while  the  white  areas  rep¬ 
resent  the  “off”  value.  From  Figure  4.7,  the  extinction 
model  predicts  the  25  km  case  to  have  essentially  no  ex¬ 
tinct  locations  in  the  plume  shear  layer.  For  30  km,  the 
model  predicts  that  much  of  the  plume  is  extinguished 
with  burning  allowed  in  the  plume  farfield.  At  35  km, 
the  plume  flame  is  extinguished.  This  is  an  interesting 
result  because  plots  of  the  local  strain  rate  show  that  it 
is  falling  with  increasing  altitude  due  to  falling  turbu¬ 
lence  levels.  This  may  be  seen  in  Figure  4.9  which  is  a 
plot  local  strain  rate  across  the  plume  at  a  streamwise 
distance  of  a :/xrefi  =  2.  The  values  of  s  in  the  plume 
shear  layer  is  seen  to  drop  dramatically  with  increasing 
altitude.  This  fact  would  at  first  glance  appear  to  rule 
out  the  strain  rate  extinction  mechanism  for  afterburn¬ 
ing  shutdown.  However,  the  reaction  rates  (or  small  scale 
Damkohler  number)  are  decreasing  with  altitude  faster 
than  the  local  strain  rate  so  that  the  extinction  limit  is 
reached  producing  the  characteristics  seen  in  Figure  4.8. 

Values  found  in  the  literature  for  the  coefficient  Cs  in 
equation  (2)  range  from  0.05758  -  2.43,  with  a  value  of  , 
2.43  being  used  by  Sloan  and  Sturgess  (1994).  These 
values  are  generally  determined  by  matching  model  re¬ 
sults  for  flame  standoff  distance  with  experimental  data 
for  low  Mach  number  turbulent  jet  flames.  Conse¬ 
quently,  these  calibration  values  likely  do  not  apply  to 
high  Mach  number  compressible  flames  encountered  in 
exhaust  plumes.  Calculations  using  Cs  —  2.43,  fol¬ 
lowing  Sloan  and  Sturgess  (1994),  predicted  extinction 
for  the  entire  plume  at  all  simulated  altitudes.  With 
Cs  =  0.05758,  the  simulations  at  30  and  35  km  were 
completely  extinguished,  while  at  25  km  the  plume  was 
mostly  extinguished.  The  value  of  Cs  =  0.03  was  cho¬ 
sen  because  it  was  the  closest  value  to  the  range  found  in 
the  literature  which  produced  a  fully  afterburning  plume 
at  25  km,  in  qualitative  agreement  with  observations 
of  real  Amine  boosters  similar  to  the  generic  configura¬ 
tion  investigated  here.  To  fully  implement  this  model 
with  GASP,  however,  would  require  a  rigorous  calibra¬ 
tion  study  to  determine  the  value  of  Cs  applicable  to 
high  speed  plume  flows.  However,  the  present  value  of 
Cs  =  0.03  was  reasonable  for  the  present  purpose  of 
demonstrating  the  viability  of  the  strain  rate  extinction 
mechanism  for  afterburning  shutdown. 

To  estimate  the  effect  of  the  extinction  model  on  the 
flow  results,  Case  B  was  rerun  with  the  chemistry  turned 
on  at  axial  locations  corresponding  to  “on”  values  plot- 
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ted  in  Figure  4.8.  Figure  4.10  compares  the  temperature 
field  for  the  extinction  model  case  with  Case  A  results  at 
30  km.  The  extinction  model  clearly  has  a  major  impact 
on  the  combustion  characteristic  within  the  plume.  The 
model  shows  a  large  delay  distance  before  a  rapid  rise  in 
temperature  occurs  as  turbulent  strain  rates  drop  below 
the  extinction  limit.  Conversely,  Case  A  results  show  a 
somewhat  gradual  rise  in  temperature  along  the  plume 
as  afterburning  is  initiated. 

The  modified  temperature  distribution  for  the  strain 
rate  extinction  model  seen  in  Figure  4.10  has  a  major 
impact  on  the  spatial  distribution  of  radiant  intensity 
as  seen  in  Figure  4.11.  Also  included  in  the  figure  are 
the  results  for  Case  A  at  30  km.  This  figure  shows  that 
the  delay  in  ignition  resulting  from  the  extinction  model 
produces  a  long  delay  in  high  levels  of  radiation  inten¬ 
sity,  significantly  changing  the  character  of  the  spatial 
emission.  Integrated  station  radiation  was  also  strongly 
affected  by  the  extinction  model  as  seen  in  Figure  4.12. 
The  extinction  model  causes  the  station  radiation  to  ex¬ 
hibit  a  sharp  rise  once  plume  ignition  occurs  and  sig¬ 
nificantly  reduces  the  peak  value  that  is  realized.  To¬ 
tal  intensity  plots  (not  shown)  also  show  the  model  to 
dramatically  increase  the  drop  off  rate  of  the  afterburn¬ 
ing  shutdown  event.  These  are  significant  findings  be¬ 
cause  the  dramatic  impact  of  the  strain  rate  extinction 
model  seen  in  the  previous  figures  will  have  a  signifi¬ 
cant  impact  on  the  total  heat  transfer  to  the  booster 
base  region  and  the  development  of  missile  defense  sys¬ 
tems.  These  results  also  demonstrate  the  viability  of  this 
previously  unrecognized  phenomena  occurring  in  rocket 
exhaust  plumes. 


5  Conclusions 

A  computational  study  was  undertaken  to  assess  the  va¬ 
lidity  of  physical  mechanisms  believed  to  contribute  to 
afterburning  cessation  in  fuel  rich  missile  exhaust  plume.- 
Several  mechanisms  were  investigated:  1)  a  relaminar- 
ization  phenomenon,  2)  a  Damkohler  number  effect  and 
3)  a  classical  flame  extinction  mechanism.  Analysis  of 
the  simulation  results  lead  to  the  following  conclusions: 

1.  The  relaminarization  mechanism  for  afterburning 
shutdown  was  found  to  be  implausible  for  the  fol¬ 
lowing  reasons.  First,  hot  plume  core  gases  and  high 
temperatures  in  the  plume  shear  layer,  due  to  vis¬ 
cous  energy  dissipation,  provide  an  ignition  source 
for  unburned  fuel  to  high  altitudes.  As  a  result, 
a  purely  laminar  plume  will  burn  at  much  higher 
altitudes  than  a  turbulent  one.  Second,  calculated 
plume  shear  layer  Reynolds  numbers  for  the  high 
altitude  case  show  the  layer  to  be  well  above  the 
transition  Reynolds  number  after  afterburning  shut¬ 
down  has  occurred. 


2.  The  Damkohler  number  mechanism  for  afterburn¬ 
ing  shutdown  was  found  to  be  plausible  and  indeed 
the  only  shutdown  mechanism  modeled  within  most 
commercially  available  codes.  This  results  from  the 
assumptions  used  to  model  the  mean  reaction  rate 
in  the  species  conservation  equations.  This  mech¬ 
anism  was  found  to  be  responsible  for  the  gradual 
drop  off  rates  of  plume  radiation  produced  by  the 
GASP  code.  The  drop  off  rate  produced  by  this 
mechanism,  however,  will  be  highly  dependent  on 
the  ignition  characteristics  of  the  fuel  species  within 
the  plume.  However,  it  seems  unlikely  that  chem¬ 
istry  effects  alone  could  produce  the  steep  drop  off 
rates  which  have  been  observed  for  many  systems. 

3.  The  classical  flame  extinction  mechanism  was  also 
found  to  be  a  plausible  afterburning  shutdown 
mechanism.  Overlay  calculations  employing  this 
mechanism  produced  dramatic  changes  in  the  spa¬ 
tial,  station  and  total  radiation  emissions.  This 
mechanism  was  also  found  to  significantly  enhance  # 
shutdown  and  offers  a  likely  explanation  of  the 
rapid  shutdown  behavior  observed  for  many  sys¬ 
tems.  This  is  a  significant  finding  because  strain  in¬ 
duced  flame  extinction  is  a  previously  unrecognized 
phenomena  occurring  in  rocket  exhaust  plumes  dur¬ 
ing  afterburning  shutdown. 
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Figure  3.1.  Schematic  of  the  axisymmetric  rocket 
geometry;  djdi  -  2,  d'/di  =  2,  d/dj  =  50. 


Figure  3.2.  Species  mole  fraction  profiles  across  the 
engine  exit  plane.  R  is  the  nozzle  exit  radius. 


Figure  3.3.  Streamwise  velocity  across  the  engine  exit 
plane.  R  is  the  nozzle  exit  radius  and  u„  is  the 
centerline  velocity. 
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Figure  3.4.  Temperature  profile  across  the  engine  exit 
plane. 


(a) 


(b) 


Figure  4.1.  Temperature  contours  for  Case  A  at  25  km,  (a)  farfield  view,  (b)  nearfield  view.  The  downstream 
direction  is  to  the  right. 


Figure  4.2.  C02  mole  fraction  contours  for  Case  A  at  25  km,  (a)  farfield  view,  (b)  nearfield  view.  The  downstream 
direction  is  to  the  right. 
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Figure  4.4.  Predicted  spatial  distribution  of  radiant  intensity  as  a  function  of  altitude  for  (a)  Case  A  ,  (b)  Case  C  and  (c)  Case  D. 
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Figure  4.6.  Reynolds  number  along  the  plume  shear  Figure  4.7.  Damkohler  number  variation  with  altitude 

layer  for  Case  A  at  35  km .  for  the  Case  A. 
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Figure  4.8.  Contour  plot  of  the  strain  extinction  model  binary  switch  for  Case  B  at  altitudes  of  (a)  25,  (b)  30  and  (c) 
35  km.  The  black  region  represents  “on”  values  and  the  white  represents  “off’  values.  The  x  coordinate  in  each  plot 
has  been  scaled  by  xrefl. 


Figure  4.9.  Predicted  small  scale  strain  rate  across  the 
plume  at  x/xrefj  =  2  for  Case  A. 
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Figure  4. 1 0.  Comparison  of  temperature  contours  at  30  km  for  Case  A  (top)  and  Case  B  using  the  strain  extinction 
overlay  technique  (bottom).  The  streamwise  axial  coordinate  has  been  scaled. 


Figure  4.1 1.  Spatial  distribution  of  radiant  intensity  for  (a)  Case  4  and  (b)  Case  B  using  the  strain  extinction 
overlay  technique. 


Figure  4.12.  Station  radiation  prediction  for  Case  A 
and  Case  B  using  the  strain  extinction  overlay 
technique. 
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